Universality, frustration and conformal invariance 
in two-dimensional random Ising magnets 
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We consider long, finite- width strips of Ising spins witli randomly distributed couplings. Frustration 
is introduced by allowing both ferro- and antiferromagnetic interactions. Free energy and spin-spin 
correlation functions are calculated by transfer-matrix methods. Numerical derivatives and finite- 
size scaling concepts allow estimates of the usual critical exponents "i/u, a/u and v to be obtained, 
whenever a second-order transition is present. Low-temperature ordering persists for suitably small 
concentrations of frustrated bonds, with a transition governed by pure-Ising exponents. Contrary to 
the unfrustrated case, subdominant terms do not fit a simple, logarithmic-enhancement form. Our 
analysis also suggests a vertical critical line at and below the Nishimori point. Approaching this 
point along either the temperature axis or the Nishimori line, one finds non-diverging specific heats. 
A percolation-like ratio 7/// is found upon analysis of the uniform susceptibility at the Nishimori 
point. Our data are also consistent with frustration inducing a breakdown of the relationship 
between correlation-length amplitude and critical exponents, predicted by conformal invariance for 
pure systems. 

v6.1 PACS numbers: 05.50.-fq, 05.70.Jk, 64.60.Fr, 75.10.Nr 



I. INTRODUCTION 

Frustration induced by quenched randomness may 
have [father complex effects on the behaviour of spin sys- 
temsEl. The Edwards- Anderson modelEI for spin glasses is 
one on which a great deal of theoretical effort has been 
concentrated, mostly on account of the basic simplicity 
of its formulation. In the version that will be of direct 
interest to us, one has Ising spins interacting through cou- 
plings of the same strength and random sign (± J Ising 
spin glass model) . For the symmetric case of equal con- 
centrations of -f and — signs, a two-dimensional [d — 2) 
system with only first-neiglibour interactions is param- 
agnetic at all temperaturesSa. Recent results, according 
to which addition of a suitable set of second-neighbour 
couplings would stabilise a low-temperature spin-glass 
phasea, have been disputed on grounds that the ordering 
thus observed is a finite-size effectB. 

However, for asymmetric distributions of ferro- and an- 
tiferromagnetic bonds one can have long-range order in 
d = 2 with only first- neighbour interactions (see Ref . R for 
references to early work) on a square lattice, for suitably 
small concentrations of frustrated plaquettesQ. Exper- 
imeutal wodc on frustrated magnets X2Cu:rCoi_a;F4, X 
= KO or Rba, has shown the existence of an ordered state 
in such conditions, though direct quantitative compari- 
son with theoretical calculations for ± J Ising systems is 
not straightforward because (1) ferromagnetic (Cu-Cu) 
bonds are weaker than antiferromagnetic (Co-Co) by a 
factor of ~ 8, and (2) the magnetism of Cu is Heisenberg- 
like with a small Ising anisotropy, which induces addi- 
tional complications such as transverse freezing. 

Estimates of the location of, and properties along—the 
critical line of the ±J model have been producedtjO. 



Further, since the pioneering work of Nishimoritj asym- 
metric Ising spin glasses have been shown to display 
quite unique features, in all space dimensionalities. 
Among these is the so-called Nishimori line (NL) along 
which the internal energy can be calculated exactly, and 
which roughly separates the regions in the temperature- 
randomness parameter space where either ferromagnetic 
or spin-glass correlations dominate. Accordingly, the in- 
tersection of the NL with the para-ferromagnetic bound- 
ary is of special significance, even in d = 2 where no spin- 
glass phase is expectejiiat non-zero temperature (T 7^ 0). 
It has been proposedMl that such intersection, to be re- 
ferred to as Nishimori point (NP), coincides with the 
multicritical point expected to exist along the boundary. 
Very recently, the NE has been investigated in d = 2 and 
4 by series analysialil, aud in d = 2 also bji a variant of 
the Chalker-CoddingtorJia network mode£j. Numerical 
values of critical exponents thus obtained are very close 
to those of percolatioiO in d = 2 (though not in d = 4EJ). 
On the other hand, d = 2 unfrustrated random Ising 
systems, such as the random-bond (i.e., bonds being J 
or rJ, < r < 1), and the diluted model, have been the 
subject of renewed interest over the past years, for two 
main reasons. Firstly, different scenarios regarding the 
universality class of these systems have been proposed: 
weak universality versus logarithmic "corrections". In the 
former, critical exponents are distinct from those in the 
pure case, but their ratios to i^, the correlatkai length ex- 
ponent, remain the same as in the pure casell3; in the lat- 
ter, the pure-system power-law critical behaviour is rein- 
forced by logarithmic divergencesEa. Secondly, the appli- 
cability of conformal invariance to random spin systems 
has ruit been put on grounds as firm as those for pure sys- 
temsEj, thus (for instance) the corresponding relationship 



between critical exponents and correlation length ampli- 
tudes needs to be checked in each case. 

Thus, a systematic study of the asymmetric ± J Ising 
model on a square lattice is of interest, not only in rela- 
tion to specific quantitative questions (such as the shape 
of the critical line and its intersection with the NP), but 
also in relation to the broader context of universality 
classes in disordered systems, singling out the effects of 
frustration. With this in mind, here we will focus on 
the following main questions: (i) to what extent, if any, 
do logarithmic corrections to pure-system behaviour de- 
scribe criticality for small degrees of frustration? (ii) does 
the connection between critical exponents and correla- 
tion length amplitudes hold in the case? (iii) |Can we 
provide evidence for (or against) the conjectureEJ that 
the critical behaviour at the d = 2 Nishimori point is 
percolation-like? These issues are addressed through the 
calculation of free energies and spin-spin correlation func- 
tions on long, finite-width strips of a square lattice. The 
rate of decay of correlation functions determines corre- 
lation lengths along the strip. We have already shown 
how averaged values of such quantities, and their nu- 
merically calculated field- and temperature derivatives, 
enable one to extijact-icritical properties of unfrustrated 
disordered modelsEiJc3. For the latter class of systems 
in particular, we gave numerical evidence in favour of 
the logarithmic corrections scenario; we also showed that 
the relationship between critical exponents and correla- 
tionJfingth amplitudes, predicted by conformal invari- 
ancecj, remains, valid provided one uses averaged corre- 
lation lengthsElra. The validity of conformal invariance 
ideas for (unfrustratftd^j-disordered g-state Potts models 
has also been verifiedEJtJ. 

This paper is organized as follows. In Sec. |l^ we outline 
numerical aspects of our calculational procedures, as ap- 
plied to the asymmetric ± J Ising model. Results for the 
phase boundary and critical behaviour above the Nishi- 
mori line are discussed in Sec. Ill, while the Nishimori 



point and the region below it are discussed separately in 
Sec. [V. Our findings are then summarized in Sec. M. 



II. CALCULATIONAL METHOD 

We have used long strips of a square lattice, of width 
4 < L < 14 sites with periodic boundary conditions 
across the strip. Only even widths were used, in or- 
der to accommodate possibly-occurring unfrustrated an- 
tiferromagnetic ground states. We compute spin-spin 
correlation functions aloagj-the "infinite" direction by 
transfer-matrix methodsBScS, extracting averaged cor- 
relation lengths. By the same methods we numerically 
obtain the free energy and its second derivatives with re- 
spect to (i) a uniform external field, which are used in 
connection with finite-size scaling (FSS) for estimating 
-f/v; and {ii) temperature, again used with FSS con- 
cepts for estimating a/iy. In order to provide samples 



that are sufficiently reppesentative of disorder, we iter- 
ated the transfer matrixH typically along 10^ (10^ near 
the NP) lattice spacings. 

At each step, the respective vertical and horizontal 
bonds between first-neighbour spins i and j were drawn 
from a probability distribution 



P(Jy) = p 5{J,j - Jo) + {l-p) 5{J,, + Jo) 



(1) 



For a square lattice the phase diagram in the T — p plane 
is invariant with respect to the symmetry p <-^- l—p; thus 
we shall restrict ourselves to 0.5 < p < 1, meaning that 
bulk antiferromagnetic order will play no part in what 
follows. 

The direct calculation of correlation functions, (ctocttj), 
goes according to Section 1.4 of Ref. |2^, with the cat- 
responding adaptations for an inhomogeneous systemEj. 
For fixed distances up to i? = 50, and for strips with 
lengths as given above, the correlation functions are av- 
erage d over a n ensemble of 10^-10^ different estimates to 
yield {(To(Jr). The average correlation length, ^"■'" (which 
carries a dependence on T, p and strip width L), is in 
turn defined by 



iaoan) ^ exp {-R/r) , 



(2) 



and is calculated from least-squares fits of straight lines 
to semi-log plots of the average correlation function as 
a function of distance, in the range 10 < i? < 50. Fi- 
nally, (,"''" is itself averaged over the different realizations 
of disordered bonds. 

In this context it must be recalled that, although in- 
sample fluctuations of correlation functions do not die 
out as strip length is increased, jOyeraged values converge 
satisfactorily^^; as done beforecll, here we make use of 
this fact to calculate error bars of related quantities. 

In contrast with-iie unfrustrated disordered models 
considered earlierOS, here the exact critical tempera- 
ture is not known as a function of p, so our first step 
was to use-iSweraged correlation lengths together with 
FSS ideasB^iSS to obtain an approximate critical curve 
Tc{p). This approach is safe because the only underly- 
ing assumption is that a second-order phase transition 
occurs, without further hypotheses on its universality 
class. In the usual phenomenological renormalisation 
recipe, used for pure systems, one looks for the fixed 
point T* of ^"''(L,T*,p)/L = ^"''{L' ,T*,p)/L' (in the 
case one would use L' = L — 2). For disordered systems 
it should be stressed that, even if logarithmic corrections 
are present in the bulk limit, the (averaged) correlation 
length at the critical point should still scale linearly with 
the strip width L, to leading orderEll. Thus, here we pro- 
duce estimates of Tc{p) by scanning a range of temper- 
atures for fixed p, and bracketing the interval for which 
^""/L appears to approach a finite value as L — > oo. The 
width of such temperature interval gives the respective 
error bar, as illustrated in Fig. |l|. Two remarks are in or- 
der in relation to this approach. Firstly, this is more con- 
venient here than the standard fixed-point search, since 



(1) intrinsic uncertainties associated to the individual ^"^ 
are amplified when estimating T*{L,L'), and (2) it is 
the extrapolation a,s L,L' ^ oo that matters in the end. 
As the range of available strip widths is not very broad, 
it is important that, for given T, the sequence of data 
from individual Us has one more point, and also slightly 
smaller error bars, than that of T*{L,L'). And, sec- 
ondly, with our procedure one already gains an insight 
into corrections to scaling: by varying the power of 1/L 
against which ^°-'" /L is plotted, one can check how bet- 
ter to produce (inclined) straight lines within the brack- 
eted temperature range singled out by the initial scan. 
It must be stressed that the location and width of the 
bracketed range itself, separating the (high-temperature) 
regime in which one is certain that ^"■'" /L -^ and that 
(low-temperature) in which ^"■'" /L diverges, are practi- 
cally insensitive to the choice of power. Indeed, though 
in Fig. H-jAfe plotted ^°-^ jL vs. L~^j^ inspired by results 
for pur^9 and unfrustrated randomo systems -, we have 
also checked that using 5,°"^ /L vs. L^^, L/^°-^ vs. L~^ 
or L~^, changes no significant digits of our extrapolated 
estimates. 

Once, for fixed p, one has an estimate of Tcij)), the 
next step is to calculate the critical free energy and its 
appropriate derivatives. This is done by evaluating the 
largest Lyapuip«-exponent A^, for strips of width L and 
length N ;» lEilH. The average free energy per site is 
//"(T) = -^Al, in units of Jq. The initial suscepti- 
bility and specific heat of a strip are then given by: 




FIG. 1. Size-scaled correlation length versus 1/L (L is the 
strip width) at different inverse temperatures, K = Jo/ksT 
in the low frustration regime: (a) p — 0.95 and K — 0.535 
(stars), 0.534 (triangles), 0.533 (crosses), and 0.532 (squares); 
(b) p = 0.92 and K = 0.637 (hexagons), 0.636 (stars), 0.635 
(triangles), 0.634 (squares), and 0.633 (crosses). Lines drawn 
through the points are guides to the eye. See Table M for 
corresponding estimates of Kc. 



Xl{T) 



d^friT) 



dh^ 



Cl{T) = 



d^ir 



(T) 



h=0 



dT^ 



(3) 



h=0 



where h is an external field coupling to the order pa- 
rameter; the size dependence of these quantities will be 
discussed below. We shall take h as uniform (ferromag- 
netic order), which is reasonable for low frustration; at 
the Nishimori point, this choice implies singling out one 
of the two scaling directions (more on this below). An 
extensive discussion of calculational details is given, for 
the specific heat, in Ref. ^. Here we recall that, since the 
derivatives are numerically obtained by calculating e.g. 
2/L(re) - MTc + ST) - hiT, - ST) with ST = IQ-^T,, 
sample-to-sample fluctuations are roughly as large as the 
difference between free energies at these three tempera- 
tures; thus one must ensure that the same configuration 
of bonds (that is, the same sequence of pseudo-random 
numbers) is used in the comparison of different temper- 
atures: free energies of the same bond geometry have 
to be subtracted. The probable errors for the free en- 
ergy differences are then much smaller than those for the 
free energies themselves. Similar procedures were used in 
a transijes-matrix study of interface energies in random 
systemst2l. The same argument applies for the suscepti- 
bilities, substituting Sh (typically of order 10^** in units 
of J) for ST. 



Finally, one should have in mind that the inverse of ^"^^ 
is, in principle, distinct from the difference between the 
two leading Lyapunov exponents, which gives the decay 
of the most probablftr-or tywcfll (as opposed to averaged) 
correlation function.L2lSt3E3a Nonetheless, for d — 2 
unfrustrated disordered Ising systeins they have turned 
out to be numerically very closeE3'LJ, at least at the crit- 
ical point (see below for remarks on low-temperature 
behaviour in the present case); significant differences 
arise only in the corrections to scaling, which are—rele- 
vant for extrapolation to the thermodynamic limita. In 
Refs. KO, the model considered here was studied with 
the aid of typical correlation lengths, ^*i'P, also calculated 
on strip geometries, but disregarding corrections to scal- 
ing; below, we will comment on some of the differences 
between our results and theirs. 

TABLE I. Inverse critical temperatures for low frustra- 
tion. The pure-system value of ^/L is A/tt = 1.2732 . . ., from 
conformal invariance. 

p Kc{p) limL^oc ^1 L 

0.99 0.4555 ± 0.0005 1.275 ± 0.015 
0.95 0.534 ±0.001 1.285 ±0.015 
0.92 0.6360 ±0.0015 1.325 ±0.015 



III. ABOVE THE NISHIMORI POINT 

We start by applying the above-described procedure 
to scale correlation lengths, for p close to 1. Figure |l| 
displays £,/L vs. 1/X^ at different temperatures, for (a) 
p = 0.95 and (b) p = 0.92. The corresponding estimates 
for the critical temperatures are shown in Table m, and 
compare rather well with those of Refs. ^ and |lj. Using 
exact and approximate data, respectively at p = 1 and 
0.99, the reduced slope of the critical curve at the pure 

= 3.25±0.11,which 



point is estimated to be ^ ,,s -j^ = 3.25±0.11, 
compares very well with the exact result, 3.20910. 
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and Tc{p), we obtain one estimate of z/^ for each disor- 
der configuration; these estimates are then averaged over 
different disorder configurations to yield the data shown 
in Fig. H, for L — 6 — 14. The trend displayed in Fig. 
is dramatically different from the one observed in the 
case of unfrustratcd disorder^ All curves (for different 
values of p) show a distinct downturn (as L increases), 
and a limiting value i/qo = 1, common to all values of p 
considered, becomes more likely. 

This should be contrasted with the case of unfrustratcd 
disorder, for which no downturn was observed, and the 
extrapolations indeed seemed to indicate a steady conver- 
gence towards a disorder dependent exponent v = J^(p) r>. 
Impure — 1; c.f. Ref. |2|. In that case, available theoriesta 
pointed towards describing those apparent exponents as 
resulting from power-law divergences with pure-system 
exponents, enhanced by multiplicative logarithmic terms; 
such expectations were latei: r.arifirnied through transfer 
matrix calculations on stripa£j'Ej. Though in the present 
case the downturn in the trend may be taken as indica- 
tive that these corrections are absent, considerable in- 
sight should be gained by trying to fit the data along 
similar lines. 

The forms of logarithmic corrections in random sys- 
tems kave been derived within a field-theoretic ap- 
proachlij, which does not explicitly account for frustra- 
tion effects. Nonetheless, inspired by our experience with 
unfrustratcd disorder, we decided to check whether m 
the present case such corrections also arise. The theoryt^ 
which successfully accounts for unfrustratcd disorder pre- 
dicts that the correlation length of the disordered Ising 
model, near the critical point, is given by 



Coo-^t-'^il + Clnil/t)]'- 



(6) 



FIG. 2. Size dependence of ul as given by Eq. (Isf), for 
p = 0.99 (triangles), 0.95 (crosses), and 0.92 (squares). Error 
bars are smaller than data points. 



We now turn to the correlation-length exponent, v. 
Since v does not appear explicitly in the expression for 
£,f"{Tc)/L, one resorts to the temperature derivative of 
the correlation, length, which can also be cast in a similar 
scaling form,Ej 



Mi 



da 



dt ^ '' L ' 



(4) 



with t = {T — Tc)/Tc, and Q is a finite-size scaling func- 
tion. Assuming a simple power-law divergence ^oo ^ t~'^ , 
- i.e., ignoring, for the time being, less-divergent terms 
such as logarithmic corrections - we obtain the estimates 
for systems of sizes L and L — 2: 



1 _ \n{p.L/HL-2)T=T, 

'^ ^ ln(i/L - 2) 



1, 



(5) 



where the derivatives are calculated numerically at the 
extrapolated (i.e., L —> oo) value of Tc- For fixed L 



for the infinite system, where v = 1, C is a disorder- 
dependent positive constant, and i> — 1/2; for C — one 
recovers pure-system behaviour. As discussed in Ref. |2l|, 
logarithmic corrections do not show up in the correlation 
length for finite systems, but in its temperature deriva- 
tives; at criticality, i.e. t{p) = 0, Eq. (^) becomes 



L2 



(1-^ InL)" 



(7) 



where A is some disorder-denfindent constant. While in 
the Dotsenko-Shalaev theoryt£l, £> was predicted to be 
1/2, here we allow it to be determined from an analysis 
of the data: it is chosen in such a way that, for fixed 

p, data for [/i^/L^] versus hiL lie on a straight line, 
for the largest system sizes. Figure shows our results 
for p — 0.95 and 0.92: we see that an attempt to fit our 
data to these logarithmic corrections would require an 
unlikely variation from i> = 0.75 for p = 0.95 to ;/ = 100 
for p = 0.92. Attempts to fit the data to other forms, 
such as powers of In L, turned out to be equally unsuc- 
cessful. The conclusion is that logarithmic enhancements 
play no role in the bulk correlation length for frustrated 



disorder (at least in a simple, clearly-defined way as in 
the unfrustrated case). 
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FIG. 3. Fits of Eq. (|) to determine v, for p = 0.95 (tri- 
angles) and 0.92 (squares). Error bars are smaller than data 
points. 
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FIG. 4. Specific heat as a function of InlnL, for p — 0.95 
(triangles) and 0.92 (squares). Error bars are smaller than 
data points. 



Turning now to the specific heat behaviour, we recall 
that in the Dotsenko-Shalaev theory, the singular part of 



the bulk specific heat per particle for the disordered Ising 
model, near the critical point, is given by 



Coo(t)^(l/Co)ln(l + Coln(lA)) 



(8) 



where Co is proportional to the strength of disorder, and 
the pure-system simple logarithmic divergence is recov- 
ered as Co -^ 0. For Co j^ and t ^ 1 a double- 
logarithmic singularity arises, whose amplitude Eq. m) 
predicts to decrease as disorder increases. For a fini|ta 
system, the usual FSS theory applied to this case yieldsE£l 



CLit = 0) ~ Ci + aln (1 + bliiL) , 



(9) 



where, similarly to Eq. (g|), 6 — > for vanishing disorder. 
We tried to check whether such forms had any relevance 
in the present case. Our results are displayed in Figure 
y, and a trend similar to unfrustrated randomness is ob- 
served: for low disorder, the specific heat increases with 
system size faster than in a double-logarithmic fashion 
(e.g., with InL); as disorder increases {p = 0.92), the 
best fit of the data crosses over to double-logarithmic 
behaviour. Though this may be interpreted as signalling 
the existence otiogarithmic corrections, such a discussion 
is rather subtleE3. At any rate, an inequivocal conclusion 
to be drawn from our data is that the specific heat di- 
verges as L — > cxo. Accordingly, this enables us to set 
a > in the hyperscaling relation dv = 2 — a, to obtain 
the condition u < 1. This condition, together with the 
absence of logarithmic corrections for ^, and the down- 
turn in the sequence of estimates for v, lead to a scenario 
of J/(p) = 1, as in the pure case. 

In order to build up a fuller picture of the low- 
frustration regime, we turned to an alternate quantity, 
the susceptibility. The ratio ^/ly can be obtained in the 
usual way. 



In [xl/xi 



-2Jt, 



ln[L/L ~ 2] 



(10) 



where Tc is understood to be the extrapolated value. We 
checked for self-consistency of critical-point locations and 
properties in the following way. Firstly, the procedure we 
used above to obtain Tc{p) from extrapolations of ^^^ can 
be repeated for ^'^^P; this yields a slightly different extrap- 
olated value, Tc(^*^P). For p = 0.95, for instance, one has 
Kc{i*yp) = 0.531 ±0.001, while Kd^'''') = 0.534 ±0.001. 
The sequence of susceptibilities calculated at these es- 
timates of Kc gives rise, through Eq. ([lO|), to the data 
shown in Figure |[ One clearly sees that the Ising value 
7/j/ — 7/4 is compatible with extrapolation of data 
calculated at Tc{£,"''"), and not with those coming from 
T'c(C*^^)- We take this to mean that: (i) from suscepti- 
bility, specific-heat and correlation-length data the most 
likely self-consistent picture is one in which the critical 
behaviour is pure-Ising for low frustration; (ii) though 
very likely <^'''" and ^'''^ will eventually scale in a sim- 
ilar way, higher-order corrections still produce sizeable 
distortions in the accessible range of strip widths; {Hi) 



further, it seems that, for not very large widths, ^"^ be- 
haves more reliably. 



1.75 



1.7 



1.65 



1 — I — I — I — I — I — I — I — I — I — I — I — I — r 








0.01 0.02 



0.03 



\/U 



FIG. 5. Ratio of critical exponents, "j /v, as a function of 
I/Z/^, for p = 0.95, from Eq. (lid), calculated at Tc determined 
through the scaling of ^"" (crosses) and of (_^^^ (squares) . 



The above analysis, together with the scaling law 
r) ~ 2 — {-^/v), predicts that, for low disorder, the ex- 
ponent describing the decay of correlations at critical- 
ity sticks to the pure system value, 77 = -j. Thus, if 
the exponent-amplitude relationship of conformal invari- 
ance remains valid in the case of frustrated disorder, we 
should have Yitcll^^ ^l[Tc{p)]/L = l/vrry = 1.273 .... We 
obtain an estimate of ViTciL-toc £,L[Tcij))\/ L by observing 
the trend followed by the sequences of points calculated 
at Kc{p) and at Kc{p) ± 0.001, to determine the central 
estimate and its error bars (see Fig. |l|); the outcome is 
shown in the last column of Table I. In spite of the arbi- 
trariness of this approach, the error bars thus obtained 
are certainly overestimated. Nonetheless, even with such 
generous allowances, the data for p = 0.92 show that the 
conformal invariance prediction is definitely not satisfied, 
since it lies way outside the range of the error bars. As 
the critical behaviour should be the same along the crit- 
ical line (at least within the low disorder region), we are 
led to conclude that, unlike the case of unfrustrated dis- 
order, the exponent-amplitude relationship of conformal 
invariance breaks down in the case of frustrated disorder. 

At p ~ 0.89, the transition vanishes abruptly, meaning 
that we do not find any temperature at which correlation 
lengths scale linearly with strip width. In Ref. |l^, it is 
found that the typical correlation lengths ^^^^ still scale 
linearly with L at suitably low T for a broader range of 
p-variation, along a line that significantly departs from 
the vertical on a p — T diagram; however, they find a 



maximum in ^*^p as a function of temperature for fi- 
nite values of T. This unexpected behaviour has indeed 
been observed in studies of ^*^p for (unfrustrated) disor- 
dered and random-field Ising systemstj. No similar peak 
structure occurs when we investigate average correlation 
lengths; instead, these vary monotonically and diverge 
only as T ^ 0, consistent with the fact that a strip is es- 
sentially one-dimensional. Though the authors of Ref. \W^ 
acknowledge that such maxima at finite T are unphysical, 
they assume that their data still are reliable above the 
peak temperatures, and interpret the corresponding part 
of their low-temperature results as marking the boundary 
between a random-antiphase state and the paramagnetic 
regime, extending as far as p ~ 0.8 . We have not found 
any evidence for this phase from our treatment. Strictly 
speaking, this means only that the expected signature 
of the corresponding second-order phase transition does 
not show up when averaged correlation lengths are con- 
sidered. At present we are unaware of why it should 
be so, and whether it means that the random-antiphase 
state is not present at all. 



1.9 



1.85 




1.7 - 



1.65 - 



T — I — I — I — I — I — I — I — 1 — I — I — I — r 




J I L 



J I L 



J I L 



0.02 0.04 0.06 
1/L2 



FIG. 6. Size-scaled averaged correlation length versus 1/L^ 
along the Nishimori line: p — 0.892 (uptriangles) , 0.891 (di- 
amonds), 0.890 (stars), and 0.889 (down-triangles). For each 
p, the corresponding temperature is given by Eq. (O). Error 
bars are smaller than data points. 



IV. AT AND BELOW THE NISHIMORI POINT 



The Nishimori line is given byt3 

exp(2Jo/T)=p/(l-p), 



(11) 



and our first task is to determine its intersect with the 
critical curve, which is done as follows. For each T, we 



extract ptriai from Eq. (g_l|), and calculate C^lPtriai); this 
procedure is repeated for different system sizes, so that 
a sequence of estimates, £,l/L^ is produced. Figure g 
shows our data thus obtained, and two different trends 
can be clearly observed: curves for p = 0.892 and 0.891 
display an upward curvature, whereas those for p = 0.890 
and 0.889 are bent downwards. Assuming a monotonic 
behaviour (as L -^ cx)) of £,/ L, any curve outside the 
interval [0.8900, 0.8910] will certainly not stabilise to a 
constant value for larger L. Our central estimate for the 
NP is therefore just the midpoint along the confidence 
interval (or, one might say, along the complementary of 
the non-confidence domain): 



PN = 0.8905 ± 0.0005, Tn = 0.954 ± 0.002, 



(12) 



where Tn follows from Eq. ( |llD . Our estimate for the lo- 
cation of the NP should be compared with those coming 
from: series work on the NLEj, giving pjv — 0.886±0.003, 
T/v = 0.975 ± 0.006; zero-temperature calculations, to- 
gether with a no-reentrance assumption, giving p^ — 
0.896 ± 0.001 |Os pw = 0.894 ± 0.002 (depending on de- 
tails of the fit)Ea; exact combinatorial worktj p^ ~ 0.885 
[error bar presumably ~ 0.005 (our estimate)]; &utX 
Monte Carlo analysis of non-equilibrium relaxation,c3 
PN = 0.8872 ± 0.0008. 
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FIG. 7. Size-scaled averaged correlation length versus 1/L^ 
below the Nishimori line: p = 0.892 (uptriangles) and 0.889 
(down-triangles) at T = 0.4. Error bars are smaller than, or 
of similar size to, data points. 



Once the Nishimori point has been accurately de- 
termined, we can check the region below for reentrant 
behaviour. Evidence has been presented very recently 
(for Potts spin-glasses on hierarchical lattices) that there 



appears to be no fundamental reason why reeatrances 
should be ruled out in thermodynamic systemsE3, thus 
this is a matter worthy of consideration. We examine the 
size dependence of ^"^ at a temperature T = 0.4 < T/v/2 
and at concentrations slightly away from pN, p ~ 0.889 
and 0.892; the results are displayed in Fig. 0. The curve 
corresponding to p > pN seems to diverge as L — > oo, 
indicating that the {T,p) point lies within the ferromag- 
netic region; the one corresponding to p < pN seems 
to vanish as L — > oo, indicating that the correspond- 
ing {T,p) point lies outside the ordered phase. Given 
that both values of p are very close to pn, we interpret 
this as an indication of absence of reentrant behaviour, 
and that the critical line Tc{p) is vertical at and be- 
low the Nishimori po iutu This is consistent with theo- 
reticai.£DtLS|iderational30 and with extensive numerical 
worklljO'Ea specifically aimed at the two-dimensional ± J 
Ising model. As a result, we assume that the scaling di- 
rections at the NP are, respectively, tangent to the crit- 
ical curve (thuSjj^iiLEfily temperature- like) and along the 
Nishimori lineB^Qo. 

In order to discuss critical exponents, we note that 
the numerical evaluation of temperature derivatives in 
Eq. (|) implies that 5T < 0.001 at the NP; since this 
is of the same order as the estimated error bars in T/v, 
we shall sit at our own central estimates, Eq. (|l2[) , and 
measure the corresponding 6T from there. 

While for scaling along the tangent {i.e., pure 
temperature-like) axis, the considerations on the need to 
subtract free energies of the same bond geometryEj are 
identical to those quoted above, a subtlety arises when 
considering variations along the Nishimori line, where a 
temperature change implies a change in p as well. From 
Eq. dn)), one has {dp/dT)p^^TM ~ 0.21. For the free 
energy calculation on what is supposed to be a given 
sample, the use of the same pseudo-random number se- 
quence for T and T ± ST with the typical ST = 0.007 (to 
be explained below) means that roughly fourteen bonds 
in 10,000 will reverse sign. We have assumed that this is 
the meaning of "using the same sample" along the Nishi- 
mori line. While in principle the bond reversals tend to 
increase in-sample fluctuations, results are manageable 
(albeit with relative error bars '^ three orders of magni- 
tude larger than those for derivatives along the pure-T 
direction), no doubt owing partly to dp/dT being small 
at the Nishimori point. We used a relatively large ST, 
compared with scaling along the pure-T axis, because 
of the need to compromise between fluctuations coming 
from the in-sample analysis described above, and the cor- 
responding actual variations of the free energy, used to 
approximate the derivative. 

We have tentatively interpreted the derivatives along 
the NL as specific-heat-like. Accordingly, we have ap- 
plied ESS to the finite-size specific heats along both scal- 
ing directions in order to find estimates of (a/iy); in both 
cases the specific heat clearly does not diverge as L ^ oo. 
Tangent to the boundary line we have found that at- 
tempts to fit our data to the form 



Cl — Co 



aL 



(?).. 



(13) 



yield a much smaller (four orders of magnitude) chi- 
square for {a/v)triai = —1-5 than for {a/v)2 — —1-1 
(the latter is extracted from v ~ 2.2 of Ref. nfl plus the 
hyperscaling relation dv = 2 — a); Figure || shows that 
our data fit neatly into a single-power form, i.e., correc- 
tions to scaling seem of little relevance in the case. Along 
the Nishimori line our data do not give a satisfactory be- 
haviour of chi-square for any sensible fitting to Eq. (13): 
varying a/v between —2 and —0.5 does not change chi- 
square significantly, and this persists even when correc- 
tions to scaling are accounted for. Thus we are not inj-a 
position to compare these data to the percolation valuetZl 
{a/v)p = —1/2. In Table || we display our raw data, so 
readers can reproduce the analysis quoted above, and try 
alternative procedures of their own devising. 
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FIG. 8. Size dependence of the specific heat at the Nishi- 
mori point, resulting from the best fit of ct/u = —1.5 in Eq. 
(ISl). The intersect with the vertical axis is at a finite value, 
0.1762. Error bars are smaller than data points. 



TABLE II. Second derivatives of free energy at the NP, for 
strip widths L = 4 — 14. T: along the temperature axis; NL: 
along the Nisihmori line (see text) . Uncertainties in last digits 
are given in parenthesis. 



L 


T 


NL 


4 


0.192183(7) 


1.558(38) 


6 


0.184854(10) 


1.704(44) 


8 


0.181824(13) 


1.772(39) 


10 


0.180245(7) 


1.846(30) 


12 


0.179288(25) 


1.846(67) 


14 


0.178667(20) 


1.859(73) 



We have thus turned to calculating the uniform sus- 
ceptibility; as it couples to a ferromagnetic order param- 
eter, the corresponding value of 'y /v is related to crit- 
icality upon crossing of the ferro-paramagnetic bound- 
ary [i.e. along the Nishimori line). In Fig. H we show 
(7/;/)^, calculated from Eq. ([lO|), with Tc = Tm, as a 
function of 1/L^. The extrapolated value 1.80 ± 0.02 
(where the estimated error bars are subjective, but cer- 
tainly conservative) compares favourably with {'y2v)p = 
43/24 = 1.7917... of percolation. Series worktfl gives 
7 = 2.37 ± 0.05 and v = 1.32 ± 0.08, which yield 
7/z/= 1.80 ±0.09. 

Finally, correlation lengths obtained from the decay of 
spin-spin correlations (which therefore couple to a fer- 
romagnetic order parameter) give an extrapolation of 
$°"/L to 1.75±0.05 at the NP (see Fig. |), rather different 
from the percolation {t^rip)^^ = 1.5279... From experi- 
ence for low frustration, as described above, we interpret 
this as signalling a breakdown of the exponent-amplitude 
relationship, rather than indicating that the transition is 
not in the percolation universality class. 
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FIG. 9. Ratio of critical exponents, 
1/Z/^, at the Nishimori point, from Eq. 
percolation value (see text). 
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V. CONCLUSIONS 

We have studied the asymmetric ± J Ising model on a 
square lattice, by means of transfer matrix calculations of 
several quantities on long, finite-width strips. Firstly, use 
has been made of a configurationally averaged correlation 
length, which is distinct from the typical (or most proba- 
ble) correlation length, used in previous transfer matrix 



studies of the same model. We have shown that an intrin- 
sically self-consistent picture can be obtained by the use 
of the former quantity, while (at least for the strip widths 
within reach) it seems that higher-order corrections to 
scaling may distort analyses based on the latter. Points 
on the critical curve Tc{p) were then obtained as those for 
which ^"■'" /L approached a finite value as L ^ cx); to the 
best of our knowledge, the estimates for Tc at p = 0.99, 
0.95, and 0.92 thus obtained are the most accurate to 
date. Secondly, the critical behaviour along the critical 
line has been discussed through the analysis of d^^^/dT, 
as well as in terms of other quantities, such as the specific 
heat and the zero-field susceptibility. 

The following picture emerged from our analysis. 
Above the Nishimori line, the correlation length and the 
susceptibility appear to diverge with power laws, with 
the same exponents as in the pure case, ly — 1 and 
7 = 7/4; logarithmic corrections (i.e., enhancements) do 
not seem to play a role in the behaviour of these quan- 
tities. We were also able to establish that the specific 
heat diverges, though at most logarithmically. Further, 
the exponent-amplitude relationship of conformal invari- 
ance breaks down as a result of frustration. These results 
are in marked contrast with the case of unfrustrated dis- 
order, for which logarithmic enhancements were needed 
in order to explain an apparent disorder dependence on 
estimates for v, and the conformal invariance prediction 
applies. 

The intersection of the Nishimori line (NL) with the 
critical curve (T^^p^) has been determined, near which 
the critical behaviour was analysed; our estimates for 
(TjVjPtv) are also the most accurate to date. Approach- 
ing this Nishimori point either along a temperature-like 
direction or along the NL, one finds non-diverging spe- 
cific heats; while for the former we were able to extract 
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(a/i/)^ 



— 1.5, for the latter we could not find reli- 



able fits. However, analysis of the uniform susceptibility, 
which probes the phase transition along the Nishimori 
line, showed percolation-like behaviour, in the sense that 
"f/iy is very close to the percolation value. Conformal in- 
variance is also absent at the Nishimori point. Further 
work is clearly necessary in order to fully elucidate all 
the subtleties related to the multicritical behaviour at 
the Nishimori point. 

Finally, we found no signature of a reentrance in the 
phase diagram below the Nishimori point; instead, the 
critical curve below this point seems to be parallel to the 
temperature direction. 
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